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Abstract 

We study modulational instability of matter-waves in Bose-Einstein condensates 
(BEC) under strong temporal nonlinearity- management. Both BEC in an opti- 
cal lattice and homogeneous BEC are considered in the framework of the Gross- 
Pitaevskii equation, averaged over rapid time modulations. For a BEC in an op- 
tical lattice, it is shown that the loop formed on a dispersion curve undergoes 
transformation due to the nonlinearity- management. A critical strength for the 
nonlinearity-management strength is obtained that changes the character of insta- 
bility of an attractive condensate. MI is shown to occur below(above) the threshold 
for the positive(negative) effective mass. The enhancement of number of atoms in 
the nonlinearity-managed gap soliton is revealed. 
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management, optical lattice, gap soliton 
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1 Introduction 



The phenomenon of modulational instability (MI) of nonlinear plane waves 
under different types of management of the system parameters has been the 
subject of intensive research over the last years pQ. Main emphasis was given 
to dispersion- management and nonlinearity-management. In nonlinear optics 
strong and rapid modulations of the fiber dispersion is achieved by periodic 
arrangement of fiber spans with alternating sign of the dispersion. Dispersion- 
managed solitons supported by such a system have essential advantages over 
conventional optical solitons for long distance communication purposes [2|3lH] . 
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Modulations of the nonlinearity is a challenging problem also in fiber ring 
lasers and in generation of Faraday waves in Bose-Einstein condensates (BEC) 
[qToITTIISIIU] . MI in the form of Faraday waves can be observed both in attrac- 
tive and repulsive condensates. Recent observation of the MI in optical media 
resulted from the periodic modulation of the nonlinearity in the evolution vari- 
able, confirms the existence of parametric resonances in the MI growth rate 
[6lfT0] . Faraday waves (parametrically excited waves) in a BEC emerging from 
temporal periodic variation of the atomic scattering length have been studied 
in [9J. Such type of modulations can be achieved by variation of the external 
magnetic field near Feshbach resonances (FR). The corresponding technique is 
known as FR management. In the Gross-Pitaevskii equation this corresponds 
to a temporal variation of the mean-field nonlinearity, i.e. to the nonlinearity- 
management. MI in a harmonically trapped BEC under FR management has 
been investigated in 



Recently the strong dispersion- management has been applied to the dynamics 
of nonlinear periodic waves, namely cnoidal waves, in optical fibers [T3|TH] . 
In these works the existence of dispersion-managed cnoidal waves and strong 
deviation of the stability borders of these waves from the ones of standard 
cnoidal wave solutions of the nonlinear Schrodinger equation (NLSE) have 
been established. Extension of the stability regions of some types of nonlinear 
periodic waves can be due to the different scenarios for the onset of MI of 
the background plane waves. Adiabatic FR management for cnoidal waves in 
optical lattices has been considered in rT51TTB] . The case of strong nonlinearity- 
management remains unexplored. 

The strong nonlinearity-management may be an effective tool for stabilization 
of matter- wave solitons in mult i- dimensional attractive BEC [17|ll8|ll9ll20|l21|l22ll23lf24 
In the context of nonlinear optics such stabilization mechanism was first dis- 
cussed in [2"Tf2"8] . The phenomenon of MI is particularly important for gener- 
ation of soliton trains in BEC with controlled spatial arrangement (repetition 
rate). MI of BEC in linear and nonlinear optical lattices in the absence of 
time-periodic nonlinearity-management has been investigated in our recent 
work [29] . Here we consider both the MI of a homogeneous BEC and MI 
of a BEC loaded in an optical lattice under FR management. The gap soli- 
ton structure existing in a BEC with the zero background scattering length 
(a s b = 0) has been investigated in Ref. [30]. The couple- mode theory can be 
used to analyze MI of nonlinear plane waves in an optical lattice subject to 
FR management. In our investigations particular interest we will paid to the 
properties of loop structures emerging in the band gaps (forbidden band). 

In the present paper we investigate nonlinear dispersion relations and the pro- 
cess of MI in a BEC under strong temporal nonlinearity management (SNM). 
The outline of the paper is as follows. The mathematical model is formulated 
in Section 2. MI in a homogeneous BEC under SNM is considered in Section 
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3. The nonlinear dispersion relation and loop structures for BEC in an optical 
lattice under SNM are analyzed in Section 4 using the coupled-mode theory. 
This section also includes the regions of MI found in different areas of the 
band structure. The properties of gap solitons are investigated in Section 5. 
Section 6 is devoted to details of our numerical procedure. In the final Section 
7 we summarize our main results. 



2 The model 



Let us consider a BEC under temporal Feshbach resonance management when 
the scattering length a s varies in time. Then an elongated BEC can be de- 
scribed by the quasi-lD GP equation with a periodic potential ( optical lattice) 
and the time- dependent management of the coefficient of nonlinearity 

h 2 

ihifh = -T^-i^xx + V(x)^ - g 1D (t)\^\ 2 ^, (1) 
2m 



where gmit) = 2ha s (t)u± is the mean field nonlinearity coefficient, uj± is the 
transverse oscillator frequency and V(x) = Vocos 2 (foc) is an optical lattice 
potential, dx\%l)\ 2 = N, N is the number of atoms. In dimensionless units 
we have 



V _, h 2 k 2 l2ha s uj ± , , , 

x -> kx, t -> u R t, e = ——, E R = — — , uj R = E R ?i, u = W — ipe . 

2E R 2m V E R 

Eq. (pQ) takes the form of the NLSE with varying in time mean field nonlinearity 
coefficient 

iu t + u xx + 7(t) \u\ 2 u — 2e cos(2x)m = 0, (2) 

where j(t) describes the strong nonlinearity-management and has the form 

/ \ 1 

7(0 = 7o + -7i ( - J , / li{r)dr = 0, r = -,//< 1. (3) 



This model has been considered in recent papers [To]ll6|3Uj . Specifically, in 
works P^lfTB] the evolution of nonlinear periodic waves under adiabatic time- 
variation of the scattering length has been studied and a possibility of gen- 
eration of a train of solitons by such a management scheme has been shown. 
Properties of gap solitons under the strong management of nonlinearity were 
analyzed based on the coupled mode system of equations in [30] • In this work 
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the gap soliton solutions and their stability for the case 70 = were investi- 
gated. Here we will study MI of nonlinear plane waves in a BEC (without and 
with an optical lattice) under SNM, as well as properties of gap solitons in the 
model ([2]) for nonzero value of 70. In particular we will analyze the possibility 
of enhancement of number of atoms in the gap soliton under SNM. 

In deriving averaged equation we follow the works [30|l3"T] and use the trans- 
formation 

1 It 

u(x,t) =e i7 - l(t)H Vx,t), 7_i(r) = / ^T r )dr' - J J ^dr'dr. (4) 

00 

Supposing the parameter \i to be small (that corresponds to high frequencies 
of modulation) unknown function v can be expanded in series as 

V = W + fJLVx + fi 2 v 2 + (5) 

where unknown w is a slowly varying function. Using transformation (j3J) and 
expansion (jSJ) in governing equation (j2J) with posterior averaging over the 
period of rapid modulation, we arrive at the following averaged equation for 

w 

i.o , . Or/iiO\ 11D 

iw t + w xx + r Yo\ w \ w — 2e cos(2x)w + a [2(\w\ ) X x\w\ + 

((\w\ 2 ) x )> = 0. (6) 

Parameter a is defined as a 2 = Jq 'y^dr. For particular case of sinusoidal 
modulations 71 = hsm(ut) we have a 2 = h 2 /(2u 2 ) ~ 0(1) (u = For 
the step-like modulation with the same amplitude h and frequency uj we have 
a 2 = h 2 /uj 2 . 

This form of averaged equation can be also obtained for the case of the weak 
nonlinearity management when 7 = 7o + li{t/n), with cr 2 C 1 [5T152] . 



3 Modulational instability of nonlinear plane wave in a homoge- 
neous media 

Now let us consider the case when the optical lattice is switched off, i.e. e = 
in Eq. (j2j). The MI of a nonlinear plane wave w = Aexp(i(7oA 2 t) can be 
explored using the linear stability analysis, i.e. looking for the solution in the 
form 

w = (A + i/j(x,t)) exp[i7 A 2 t], ip < A. (7) 
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We have the following equation for ip 

+ 1>xx + 7o^ 2 (^ + r) + ^ 2 A\^ XX + i)* xx ) = 0. (8) 

Representing if> = %p r +i%pi and performing Fourier transformation ip r (ipi)(x, t) = 
J dku(v)(k,t) exp(ikx) we get the dispersion relation 

p 2 = k 2 [2 lo A 2 - (1 + Aa 2 A 4 )k 2 ]. (9) 

Instability region corresponds to the condition p 2 > 0. Thus we obtain 

k 2 < 27 ° A * ■ (10) 
~ 1 + Aa 2 A 4 v ; 

The maximum of the MI gain is achieved at the value of the wave number 



Maximal value of the MI growth rate is 

Pc = 7TTW' (12) 



Thus we find that under the temporal nonlinearity management the MI growth 
rate is decreased by a factor of vl + Aa 2 A 4 . Such decrease of the gain is 
due to the defocusing effect induced by the nonlinearity management. This 
observation explains the stabilizing role of the strong nonlinearity management 
in a higher dimensional attractive BEC [17,18,33 f3i] . 

Numerical simulations of the ID GP equation (j2J) with a strong nonlinearity 
management confirm these predictions. In Fig. [1] we plot the MI gain versus 
the wave number of modulations k for three different cases with 70 = 1 and 
<jj — 10: (a) when the nonlinearity- management is absent, a 2 = and when 
the management is present (b) a 2 = 0.125(A = 5), (c) a 2 = 0.5 (A = 10). One 
can observe a good agreement between the theory and numerical simulations 
for the value and the position of the MI gain maximum given by Eqs. (11 II) 
and f|T2|) . In Fig. [2] we plot the profiles of the field module |w(x)| in the region 
of stability. Fig. [3] depicts the case of breakdown of the stability caused by 
increasing the strength of the nonlinearity-management, a 2 . One can see that 
modulation in an initial plane wave evolves into a train of solitons when the 
wave number of the modulation is in the region of instability. As can be seen 
from Fig. [3]^a) even moderate nonlinearity management (a 2 = 0.125) causes 
notable decreasing in the amplitude of solitons. 



5 




Fig. 1. MI gain p versus the wave number modulations k. Three curves correspond 
to the cases when: (a) nonlinear management is turned off, a 2 = 0; (b), (c) nonlinear 
management is turned on, with A = 1.2, a 2 = 0.125 and a 2 = 0.5. Filled squares 
correspond to gains obtained from full PDE simulations. 

t = 5.02 
t = 4.40 

t = 3.77 
t = 3.14 
t = 2.51 
t=1.: 
t = 1.25 
t = 0.628 



Fig. 2. Evolution of the small spatially periodic perturbation when p and the wave 
number k are in the region of stability and inequality (|10p is not fulfilled. The case 
with a 2 = 0.02, k = 2, A = 1.2 is presented. 

4 MI in a BEC loaded in an optical lattice and nonlinearity-management 



4-1 Nonlinear dispersion relation. The loop structure. 



The analysis performed in the previous section was relevant to a BEC without 
optical lattice potential. In the presence of an optical lattice the band struc- 
ture strongly affects the process of MI [35|36] . Equations of the coupled-mode 
theory for the GP equation (jH]) with shallow optical lattice have been obtained 
in [SU] • The wave function can be represented in the form of superposition of 
backward and forward propagating waves 

w{x, t) = y/e (A{X, T)e lx + B(X, T)e~ ix ) e~ il . (13) 
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time, t 



Fig. 3. Development of small spatially periodic perturbations into a soliton train 
when parameters are in the region of instability. Plot (a) depicts the field profiles 
at different times. Dashed line stands for initial small modulations at t = 0, 
dotted (solid) line is for the case, when SNM is turned off (on) at t = 20.1. Plot (b) 
depicts time evolution of the maximal value of |tt(sc)|. Dotted (solid) line correspond 
to turned off (on) SNM. The parameters are k = 0.5, A = 0.5, a 2 = (a 2 = 0.125). 

where X = ex, T = et are slow variables. Substituting this into the averaged 
equation we get the following coupled mode system of equations 



iA T + 2iA x = B- 7o(|A| 2 + 2\B\ 2 )A + 8ea 2 {2\A\ 2 + \B\ 2 )\B\ 2 A, (14) 
iB T - 2iB x = A - 7o(2|A| 2 + \B\ 2 )B + %ea 2 (\A\ 2 + 2\B\ 2 )\A\ 2 B. (15) 

In derivation of this system, the derivatives of the nonlinear terms have been 
neglected as the terms of the next order of smallness with respect to e. The 
group velocity varies in the interval — 2 < v < 2, in physical units that cor- 
responds to — vr < v < vr, vr = hk/m. This system describes two counter 
propagating waves, with the cubic self phase modulation term and cubic and 
quintic cross-phase modulation terms. The quintic cross modulation term de- 
scribes effect of the Feshbach resonance management. Note that this system 
has a similarity with the one previously considered for description of MI in the 
cubic-quintic NLSE with the Bragg grating [37] . However, as distinct from that 
model, no self-phase modulation quintic terms like |t4| 4 A and \B\ A B present 
in our model. The absence of these terms changes significantly the MI process 
in NM systems in comparison with the standart cubic-quintic NLS model. 

The plane wave solutions of Eqs. ( iT4l and ( JTBT) are looked for in the form 

A - a J(QX-QT) R _ a f i(QX-QT) 

where a=|y4| 2 + |5| 2 . The parameter / defines the weight of the forward and 
backward propagating waves. The case |/| > 1 corresponds to the domination 
of the backward wave. Substituting these expressions into the system ( TT41) and 
(1151) . we obtain nonlinear dispersion relation 
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Q = 



37o 2 
-or 



2 

(8ea 2 a 2 



1 1 + J 2 
2 / 
7o)« 2 1 



+ 



4ecx 2 a 4 



1 + / 



a + / 2 ) 2 

/ 2 , ll-/ 2 
4 



(/ 4 + 4/ 2 + X) 



./ 



(16) 
(17) 



The parameter / determines the position on the dispersion relation in Q, Q 
plane. Inspecting the dispersion relation at small a 2 one can observe that 
/ > corresponds to the upper dispersion curve and / < to the lower one. 
The velocity inside the grating is v — 2(1 — f 2 ) /(l + f 2 ) and equals to zero at 
the edges of the gap / = ±1. 

From Eqs. ([TBI) and (fT7|) one can again see a defocusing role of the strong 
nonlinearity management. We find that the effect of nonlinearity is cancelled 
if |/| = 1 and the density of BEC reaches a threshold value 



To 



4ecx 



2 ' 



Suppression of the mean- field nonlinearity in the lattice leads to enhancement 
of such an effect as tunnelling between sites. The SNM also introduces changes 
in the dispersion curves. Indeed, it is well known that the focusing Kerr non- 
linearity (attractive BEC) is responsible for appearance of a loop beyond the 
critical power [38ll39f40j on the upper curve. 

Effective nonlinear dispersion induced by the nonlinearity-management (the 
last term in Eq. (J6j) will increase the critical power necessary for appearance 
of the loop. To find this value of critical power let us consider the value of f c 
at which Q becomes zero (|/| ^ 1). We obtain that 



fc 



« 2 (7o — 8ea 2 a 2 ) 



± 



\ 



a2 {lo ~ 8ea 2 a 2 ) 



1. 



Let us consider the case of upper curve with / > and an attractive conden- 
sate 70 > 0. Then a loop appears on the dispersion curve if the power (BEC 
density) 



a\ < o? < a\, a\, 2 = ^ L ± ^ _ ] . (19) 



When a 2 = we have a well known result for the critical power [3H] ct 2 = 2/70 . 
Fig. H] depicts two branches of the dispersion relations ffTB]) and ffT7|) . The 
branches in the Q — Q plane are defined parametrically by Eq. (fl6|) and (I17p . 
Thus, each value of / defines a point in this plane. Two ranges of values / > 
and / < define upper and lower curves correspondingly. In our case the loop 
structure appears on the upper branch when a is greater than the threshold 
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" i • 1 • 1 • 1 • 1 • 1 • 1 • 1 • 1 

-2,0 -1,5 -1,0 -0,5 0,0 0,5 1,0 1,5 2,0 

Q 

Fig. 4. Loop structure in dispersion relations when a 2 > a 2 : solid (dotted) 
line and full squares (circles) are for the case when nonlinearity management is 
turned off (on). Scatter points (squares and circles) represent data obtained from 
numerical simulations. Parameters are: a 2 = and a 2 = 0.055(/i = 4) with 
a = 2.5, e = 0.08, u = 12. 

value a 2 . One can see from this figure that the loop decreases with increasing 
of the strength of the management. It should be noted that at the same time 
the band width (a distance between upper and lower branches) at Q = does 
not change. 

In the case of the defocusing Kerr nonlinearity 70 < (repulsive BEC) one 
could expect formation of the loop on the lower branch of the dispersion curve 
/ < 0. But from the condition (fl~9l) it follows that a 2 < 0. So in this case NM 
fully suppresses the loop formation. 

It is also of interest to investigate the loop structure in the case 70 = 0. This 
configuration can be realized employing the Feshbach resonance technique. It 
corresponds to the case of a BEC with the effective repulsive nonlinearity in an 
optical lattice. The loop will be formed on the lower branch of the dispersive 
curve / < when the BEC density excesses the value 



Let us discuss the physical consequences. Existence of a loop at the edge of 
the Brillouin zone reflects the superfluid character of the BEC, since we have 
nonzero velocity in the Bragg reflection condition |39)40j . It should be noted 
that in a linear system of free atoms the Bloch wave at the zone edge has 
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zero velocity. From this point of view a critical value of the SNM strength 
exists which destroys the superfluid property of the BEC in an optical lattice. 
Another possible effect is the existence of breakdown of Bloch oscillations due 
to the tunnelling into the upper band (Landau-Zeener tunnelling). The SNM 
is expected to suppress this breakdown. 



4-2 Modulational instability 



To investigate MI of matter waves in an optical lattice under SNM, perturbed 
plane wave solutions are taken in the form 



A= l^^ + 5A(X,T))eW x - nT \ (20) 

B = + SB(X, T)) e^ x -* T \ 

where 5 A and SB are unknown small perturbations of CW solutions. Substitut- 
ing these expressions into Eqs. (TH| and ffl5|) and using a linear approximation, 
we get the system of equations for 5 A and 5B 



i8A T + 2iSA x + fSA -6B + ——^[(70 - : , f2 J )(5A + 5A*) + 

2/(70 - 8ea 2 a 2 )(5B + 5B*)} = 0, (21) 

1 a 2 f 2 16ea 2 a 2 

i5B T - 2t6B x + j5B -5A+ ^^[(70 - -j— j^){8B + 5B*) + 

j( l0 -8ea 2 a 2 )(5A + 5A*)]=0. (22) 

For / = ±1 the system coincides with the one considered by de Sterke [35J with 
renormalized nonlinearity coefficient 7 r = 70 — 8ea 2 a 2 . One can see that the 
NM plays essential role in the MI process. When the nonlinearity management 
is turned off, for the case of attractive condensate the CW wave is unstable if 
the parameters follow the upper branch of the dispersion curve. On the lower 
branch the attractive BEC is modulationally stable. The repulsive condensate 
is modulationally unstable on the lower branch and stable on the upper branch. 

In the case of nonlinearity-management there exists a critical value of the 
management strength a 2 , namely a 2 = 7 /8ea 2 . If a 2 > a 2 , then the attractive 
condensate behaves as the repulsive and the modulational instability regions 
should correspond to the above described picture. 
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Looking for solutions of Eq. (I2TT) in the form 



8A(B) = C(D) cos(gX - uT) + iE(F) sin(gX - uT) 
we find the dispersion relation of the form 



(uj 2 - Aq 2 ) 2 - 2(1 - N)(uj 2 - 4g 2 ) - 1(1 +P)(u- 2q) 2 - 

/(/ + M)(^ + 2g) 2 = 0, (23) 



where 



2a 2 . 16ea 2 a 2 f 2 , Ar 4fa 2 . , 2n 

2a 2 f 2 _16eoV 



Analytical results can be obtained for the particular case |/| = 1, correspond- 
ing to the edges of the gap. We come to the equation for the frequency u 



u; 2 



V + 2 - G ± ^16q 2 (l + G) + (2 - G) , (24) 

G = G/f, f = ±1. 



Evidently, this equation coincides with the one obtained in Ref. [35] where the 
parameter G is renormalized as G = (70 — 8ecr 2 a 2 )a 2 . 

Let us analyze the condition of MI for different sets of parameters. 

1. The top of the band gap / = 1, cr 2 < a 2 (G > 0). The wave is unstable 



if the wavenumber of modulations is in the interval —y3G/2 < q < y3G/2. 
The maximal MI gain occurs at the wavenumber 



KS^r (25) 



Results of numerical simulations of the Gross-Pitaevskii equation for evolu- 
tion of the nonlinear plane wave modulations is shown in Fig. [5j The emergence 
of a train of gap solitons is observed. The reduction of the MI gain when the 
SNM is applied can be noted. 
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2. The bottom of the band gap / = — 1, a 2 < a 2 (G < 0). The condensate 
becomes unstable if G > 1 and the wavenumber satisfies the inequality 



191 > — Vg=t (26) 



3. In the case a 2 > a 2 , an attractive condensate behaves like the repulsive 
condensate under the strong nonlinearity management. We can expect mod- 
ulational instability in the case of / = —1, corresponding to the negative 
effective mass. In this case the condensate is unstable in the region of modu- 
lations with the wave numbers q 2 < 3\G\/2. 

Let us consider separately the case 70 = 0, G = — 8ea 2 a 4 . As it was shown 
in [30], near the upper edge of the gap, the gap soliton is the solution of the 
focusing quintic NLSE, while near the bottom of the gap it is a solution of 
the defocusing quintic NLSE. Fig. depicts the formation of a gap soliton 
train under strong nonlinearity management. It should be noted that when 
70 = <y 2 = 0, the soliton does not form. 



For f — —1 the instability region is |g| < y 12e<r 2 a 2 /2. For |/| ^ 1 we can 
perform analytical consideration for the case of vanishing wave numbers of 
modulations q — 0. Then in the ordinary optical lattice the gain of MI turns 
to be finite and for the MI in the normal dispersion region there exists a 
threshold in the power. In the case of the action of a SNM we find from 
Eq. (123]) that the instability occurs if 

(1 + / 2 ) 2 _ 4/( 7 o - 8eaV) 

p 1 + /2 < U - W 



For example, if 70 > and / > the MI is possible only if a 2 < a 2 < a 2 , 
where 



a 



To 



1,2 



16ea 5 



1 ± 



\ 



1 - 



8(1 + f 2 fea 2 



The MI interval on a 2 for / < can be obtained analogously. 



5 Gap soliton 



Following the works [29f30] let us study the properties of a gap soliton. The so- 
lution is sought in the form A = a(X) exp(-iClT) , B = b(X) exp(-iQT) , a = 

b*, a = JQ(X) exp(-z0pf)/2). The set of equations for Q(X), 9{X) is 
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Fig. 5. Evolution of small spatial periodic perturbations when the parameters are 
in the region of instability with 70 = 1 and / = 1 (upper branch of the dispersion 
relations (fT6|) and (fT7|) ), Plot (a) depicts the field profiles \u(x)\ at different times 
when the nonlinearity management is turned off and a 2 = 0. Plot (b) depicts the 
case when the nonlinearity management is turned on and a 2 = 0.125(/i = 5). Other 
parameters are a = 0.8, Q = 0, q = 0.5, uj = 10. Initial amplitude of modulations 
is taken to be 0.05. 
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-Q + cos(0) 



Q x = Qsm(6), 
3j Q + 2Aea 2 Q 2 . 



(28) 
(29) 



The first integral of this set is 



E = -ClQ + Q cos(#) - - 7o Q 2 + 8ea 2 Q 3 . 

Inside the gap — 1 < fl < 1. The solution for 70 7^ is difficult to be derived 
in an explicit form. What we can calculate is the peak value of gap soliton 
amplitude, the quantity, which is of interest for the experiment. For the soliton 
peak the condition Q x = is valid. Taking into account that for bright soliton 
solution E = 0, we obtain the following equation for the peak value of the 
soliton amplitude 

±1 = n + ^7oQ - Sea 2 Q 2 , 

where the signs ± correspond to 9 = and 6 = rr respectively. Peak values, 
corresponding to the bright soliton solutions, are 



Q 



32ea 



r7o 



\ 



1 + 



128e(x 2 (^T 1) 



(30) 
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t= 13.1 




t = 6.28 




t = 




Fig. 6. Evolution of small spatial periodic perturbations when parameters are in 
the region of instability with 70 = and / = — 1 (lower branch of the dispersion 
relations (JT6J) and (fT7|)). The plot depicts the field profiles \u(x)\ at different times 
when the strength of the nonlinearity-management is a 2 = 0.08(/i = 4). Other 
parameters are a = 1.5, Q = 0, q = 0.5, uj = 10. Initial amplitude of modulations 
is 0.05. 



It should be noted that when 70 = we get Q = y T l)/8ecr 2 , that coincides 
with the value obtained in |30j. Existence of two families of gap solitons has 
similarity with the ones observed in the cubic-quintic NLSE with a periodic 
potential [41]. From ( |30l) we obtain the restriction 



a 2 < 



128e(l - Q)' 



For the estimations of the experiment with e = 0.2, Q = 0.6, we obtain the 
restriction a 2 < 0.6. The defocusing role of the nonlinearity management leads 
to the possibility of increasing the number of atoms in the bright gap soliton 
in comparison with a standard gap soliton. The low nonlinearity requires the 
larger number of atoms to support soliton solution. Taking a 2 = in the 
low amplitude solution, we obtain for the peak amplitude the value Q = 
2(1— f2)/3, that reproduce the standard result for a gap soliton |42j. Expanding 
the solution fl30l) in series, we obtain 



Q ~Qo + 



64e(l -fi)V 
27t| 



The number of atoms in the gap soliton is enhanced and the enhancement 
factor is proportional to the nonlinearity map strength a 2 . For typical values 



of parameters V = O.QEr (e = 0.3), h = 3.16a;, / = 33, 



lOa^ (a 2 
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h 2 /2u 2 = 5), Q = — 1, 70 = 1 we obtain double enhancement in the number of 
atoms. It means that for the experiment with 87 Rb [13] the number of atoms 
in a gap soliton (N ~ 600) can be increased by the nonlinearity management 
up to iV ~ 1200. The increasing of number of atoms in the discrete breather of 
discrete nonlinear Schrodinger equation under weak nonlinearity management 
has been observed in numerical simulations [32] . 



6 Numerical simulations 



In numerical simulations we proceed from the governing Gross-Pitaevskii equa- 
tion (T5]). The problem is discretized in a standard way with the time step At 
and spatial step Ax so that terms Uj approximate u(jAx, kAt). More specif- 
ically, in the approximation of Eq. we have used the following implicit 
Crank-Nicholson-type scheme of second order accuracy in space and first or- 
der accuracy in time 



i{u) +1 - u)) 



At 2Ax 2 



k+l o„,fc+l i „.fc+l 

'i-i 



+ 



ecos(2x,)(u* + up 1 ) - ^T(^)I4| 2 K + up 1 ), (31) 

where the strong nonlinearity management factor j(t) is defined by Eq. (j2J), 
Xk = jAx and t^ = kAt. In our calculations the second term in Eq. (j2J) is 
chosen as 71 = hsm(ut). For this case a 2 = h 2 /(2a; 2 ). 

Since our problem deals with nonlinear plane waves, periodic boundary con- 
ditions are imposed on the governing Eq. ([2]). Eq. (l3Tj) together with the 
boundary condition up 1 = up 1 form a quasi tridiagonal set of equations for 
unknown up 1 , [j = 0, 1, 2...N] in a lattice of N + 1 points. The length of the 
lattice L is determined by the period of the periodic potential and value of 
the wave number for which the solution is sought. The set of these algebraic 
equations is solved by the modified vectorial sweep method. In actual calcu- 
lations the typical space step Ax ranged from 0.01 to 0.005 and time step At 
from 0.005 to 0.001. 

In calculations, the initial wave packet is constructed in the following way. At 
first slow component of the solution w(x,t = 0) is taken in the form of Eq. (j7j) 
or Eq. ([TBI with Eq. (|2"U1) . depending on the problem we consider. Then leaving 
only first term in Eq. (JSJ) and making use of transformation Eq. we obtain 
actual initial wave function u(x,t = 0) used in computations. 

In simulation of the loop structure (see Fig. H]) and constructing initial wave 
function, the position on the loop for given value of the wave number Q is 
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determined by choosing necessary value of the parameter /, which, in turn, is 
determined from the dispersion relation Eq. (fl6|) . 



7 Conclusion 

We have investigated the modulational instability and gap soliton formation 
in the media with Kerr nonlinearity and periodic potential. Such systems 
appear in the nonlinear optical media with Bragg grating and Bose-Einstein 
condensates in optical lattices under time-dependent Feshbach resonance man- 
agement. We considered the case of strong management and showed that in 
the case of homogeneous Kerr media under NM the gain of MI is strongly sup- 
pressed, that explains the defocusing role of the NM and thus the stabilization 
of 2D and 3D attractive BEC by this method. We have studied the nonlinear 
dispersion relation in the case of NM and showed that the loop structure is 
essentially modified by the NM. The critical value of the strength of the NM 
is shown to exist in the MI regions. In the case of attractive condensate it 
means that above the threshold an attractive BEC behaves as repulsive. The 
NM leads to a new effect of enhancement of the number of atoms in the bright 
gap soliton. The enhancement factor is proportional to the strength of the 
management a 2 . We confirmed the predictions based on the analysis of the 
averaged GP equation by direct numerical simulations of the ID GP equation. 
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